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'\^^! Abstract 

■ Comparison of solar-neutrino signals in SNO Q and Super-Kamiokande (SK) 

oo : 

^ ■ 1^ detectors results in discovery of Ue i^^i,t oscillations at level 3.1 — 3.3 a 

^ ■ This comparison involves the assumption of neutrino spectrum and a choice 

\D ' 

■ for the thresholds of detection in both experiments. In this Letter we obtain 
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an exact formula for the comparison of the signals which is valid for arbitrary 
spectra and thresholds. We find that the no-oscillation hypothesis is excluded 
p • at 3.3 cj. If the energy-dependent component of the survival probability for 



Oh! 



electron neutrinos is small as compared with the average value, i.e. in the 
case of small distortion of the observed spectrum, the oscillation hypothesis 



^ ■ can also be tested to similar accuracy. The oscillation to sterile neutrino only, 

is excluded at 3.3cr level, and oscillation to active neutrinos is confirmed at 
the same level, though with some reservations. 
Comparison of Charged Current (CC) solar-neutrino signal in SNO with Elastic Scat- 
tering (ES) signal in Super-Kamiokande (SK) [Q has revealed the presence of signal from i/^ 
neutrinos in SK detector, which evidences for neutrino oscillation. At present this proof of os- 
cillation exists at level 3.1—3.3 cr, being limited mostly by systematic errors and uncertainties 
in calculations of CC cross-section: Ue + d ^ p+p + e^ (for the recent discussion of the latter 
and the calculations of radiative corrections to the cross-section see Ref. M). The basic idea 
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of extracting the signal from another active neutrino component Ua = i^f^^r comes from the 
fact that in CC scattering in SNO experiment the flux of z/g neutrinos is measured, while the 
signal in SK is provided by all active neutrinos including u^. SNO gives the total flux of 
electron neutrinos, assuming the standard spectrum [|] as (psNO — 1-75±0.148 x 10^ cm~^s~^. 
Here and everywhere below we conservatively use the upper systematic errors and sum dif- 
ferent errors quadratically. The total neutrino flux detected by SK {also assuming the 
standard flux) is = 2.32 ± 0.085 x 10^ cm^^s^^. If one assumes that only electron and 
sterile neutrinos are present, the flux above is the flux of z/g neutrinos. If to allow the presence 
of other active neutrinos z/q, the flux above is (pi^^ + 0.1540a, where 0.154 is the ratio of cross- 
sections UaC and UeC. In any case the 3.3 cr excess (psK~^SNO — 0.57±0.17x 10^ cm~^s~^ 
signals about presence of another component of active z/^ neutrinos. Note, that both fluxes 
are obtained assuming the standard (SSM) flux and using for their calculations the different 
thresholds for electron detection. Adjusting the thresholds, the SNO collaboration arrives 
at 3.1 a flux difference 0.53 ± 0.17 x 10^ cm'^s^^ |l|]. 

An interesting and effective method of comparison of the signals from SNO and SK was 
suggested in Refs.( 0). It is based on the demonstrated property that response functions for 
SNO and SK are approximately equal at the appropriate choice of SNO and SK thresholds. 

This method has been further reflned in Ref. 0. Choosing the appropriate thresholds 
the authors verifled the presence of ocillation to active neutrinos at level 3.1 a. Among other 
interesting results the authors have obtained the value of fs, which characterises possible 
deviation of boron neutrino flux from the SSM prediction, equal to 1.03to;58. 

In Ref. it was demonstrated that in the framework of the Unifled Approach |p the 
above results correspond to the probability of z/g —>■ v^ ^ to be larger than zero at 99.89% 
CL. 

Including the results of other solar neutrino experiments and astrophysical information 
about neutrino production increase considerably the signiflcance of discovery of neutrino 
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oscillation 

In this Letter we shall derive an exact and simple formula valid for arbitrary spectra 
and thresholds. This formula will be obtained for an arbitrary process of z/g disappearance 
(e.g. oscillation accompanying by the decay of a mass-eigenstate ^^^), characterised by 
survival probability Pee{E). The oscillation, as most interesting and realistic case, will be 
analysed in detail later and anticipating it we shall keep the subscript "osc" (for oscillation) 
as a notation from the beginning. 

The basic quantity that we shall analyse is the ratio of the electron rate in SK, R^^ , 
produced by active neutrinos Ua = ^ij.,t and the total electron rate RtJ^: 

rose = <^/i?fof . (1) 

Let us introduce the following definitions. 

1. The exact flux of u^. neutrinos reaching a detector is defined as 

^uAEu) = ^B(pSSMiE,)PeeiE,), (2) 

where $b is the total boron neutrino flux at production, the value of which is not specified, 
and 4>ssm{E,^) is the spectrum at production, normalised as 

J <PssM{E,)dE, = I. (3) 

2. According to the definition of SK collaboration, the detected neutrino flux, 0f^, 
determines the total rate in SK above threshold T^f^, if the neutrino spectrum in detector is 
standard: 

^fof = 0fl/,, dT dT'RsKiTX) I dE,<PssM{E.WUEu,r), (4) 

where T' and T are the real and the measured kinetic energies of the electron, respectively, 
and Rsk(T,T') is energy resolution function of the Super-Kamiokande detector [|lT|. 
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can be interpreted as the flux of z/g neutrinos in case of v^. — > conversion, and as = 
(^y^ + 0.1540,^^ in case of conversion, wliere 0.154 is a ratio cr(z/ae)/cr(z/ee). 

3. The flux of electron neutrinos from decays, (psNO^ determined by SNO coUabo- 

■^SNO 



ration from the CC-reaction rate above threshold T^h > assuming the standard spectrum: 

rp TpmaXfrpf\ 

R'iNO = & / dT / dT'RsNoiT, T') / dE,<PssM{E,)acc{E,, T'), (5) 

Jj-SNO J JE"^^"(T') 



where RsNoiT,T') is energy resolution function of SNO 

4. On the other hand the same can be expressed via the exact flux ^^^{Ey 

as given by Eq.(0): 



RsNO = / dT / dT'RsNoiTX) / . dE,<P 

5. Finally, we define two quantities, Jj and Ji{Pee} 

MPee}= r™'rfT / dT'R,{T,T') f^^ dE,<f)ssM{E,.)Pee{EM{E.,T'), (7) 
J, = J,{P,e = 1}, (8) 

for i = SNO, SK, = cxcc for SNO and (Tj = (Tj^^e for SK. The ratio Ji{Pee}/ Ji determines 
the average survival probability {Pee)i for i =SNO, SK. This ratio is Pee in case of constant 
survival probability. 

We are now ready to obtain the exact expression for the ratio rose from Eq.(|l]). Rear- 
ranging it as rose = 1 — Ru!^ /P-tat ^ using for SK rate induced by z/g R^^ = ^sJsKiPee} 
we obtain 

"^BJsKiPee} 



'osc — J- J , yiJ) 

(PSK^SK 

where $b is given by Eqs.(|) and @ as 

= 0g^oy^^- (10) 
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We can decompose Pee{E) in the energy independent part Pee, e.g. the average value in 
one of the experiments, and a small (according to experimental data) energy dependent part 
SP{E): Pee{E) = Pee + 6P{E). Then we can re-write Eq.(|) as 

^ , <l>^S^O l + JsK{'-t}/JsK <^^s^ol + mSK 

where the angle brackets mean energy averaging. The formula (^^) is exact. It is valid in 
particular even when 6P{E) is not small in comparison with Pee- K is derived for any process 
of z/g disappearance, including the oscillation. 
The correction factor in Eg. (]TT|) , 

1 + {^}SN0 

is close to 1 in each of two cases: when (i) 5 Pee/ Pee is small or when (ii) average values 
{6Pee/Pee) for SNO and SK are close to each other. As we will see in the physically relevant 
cases both conditions work simultaneously. 

The first condition is satisfied because the observed distortion of spectrum is small. The 
correction factor (|12D is invariant relative to decomposition Pee{E) = Pee + 6Pee{E) with 
arbitrary Pge- Since the distortion of spectrum is small we always can choose the value of 
Pee such that 6Pee{E) chaugcs its sign in the middle of the observed energy interval. It 
provides the smallness of {6Pee/Pee)i- The condition (ii) will further diminish 1 — Kcon- 

In general case Kcorr can be obtained directly from observational data using 6Pee{E) for 
each energy bin. 

We shall consider now the application of Eq.(pTD for neutrino oscillations assuming that 
oscillation is the only way of Ve disappearance. 

To test the different hypotheses we shall treat rose, given by Eq.(0), as experimental 
value to be compared with a theoretical value rose = 0. The latter is the predicted value for 
the hypothesis of oscillation to sterile neutrino only, as in the case of two neutrino mixing 
(hereafter we shall refer to this case as sterile-neutrino oscillation) . In the case of oscillation 



to active-neutrinos z/^, to be referred to as active-neutrino oscillation, rose = corresponds 
to absence of — > z/^ oscillation. 

For the case Pee{E) = const we obtain 

rose = 1 - -g^ = 0.246 ± 0.0694. (13) 

rsK 

In particular, this value of the ratio is exact for no-oscillation case Pee = 1? and formally it 
is 3.54 cr away from no-oscillation value rose = 0. 

However, evaluation of the error for the ratio (I^sno/^sk does not correspond to to the 
usual connection with the confidence level, because the distribution is not Gaussian. More 
correct evaluation of the error is as follows. 

In the plane {(ps^o y 4'sk) plots experimentally measured point and la, 2a and 3 a 
contours around it (see Fig. 1). One can see that the contour 3.33 a touches the hne rose = 0. 

In the general case to obtain the value rose from Eq. (|TT|) we evaluate Kcorr in the following 
way. 

First we assume, and later verify, that correction factor for a tested hypothesis is given 
by value Kcorr ~ 1 within very narrow interval Af^corr, so that we can treat Kcorr just as a 
number. Then using 4'sno/4'sk — 0.754 ± 0.0694 we find formally for this hypothesis 

rose = 1 - 0.754 i^eorr ± 0.0694ireorr (14) 

To obtain the range of allowed values of -ft'corr one must compute the correction factors 
for all points in the parameter oscillation space allowed by all available data e.g. at la, with 
the Super-Kamiokande spectra data being the most important ingredient. 

If the starting assumption about narrow interval AKcow is confirmed, one can proceed 
further choosing the maximum value of -R'eorr from obtained range and inserting it to Eq.([T^). 
-^cor? gives the smallest value of rose with the largest error, which we conservatively take as 
the final result. 
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FIG. 1. SNO and Super-Kamiokande fluxes allowed at 1 cr, 2 a and 3.33 a levels in comparison 
with no-oscillation value rose = for two values of correction factor Kcorr = 1 and Kcorr = 1-01. 
The experimental point is 3.3 a away from rose = which corresponds to absence of active neutrino 
oscillation. 
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In both cases, sterile-neutrino oscillation and active-neutrino oscillation, a statistically 
significant deflection from rose = is a proof of the hypothesis (absence of sterile-neutrino 
oscillation, and presence of active- neutrino oscillation), though for sterile neutrinos this 
criterion is formally more strict, because rose = is the theoretical prediction of sterile- 
neutrino hypothesis. 

For sterile neutrinos -ft'eorr has been calculated [] for all points in the sterile oscillation 
parameter space (Am^,tan^^) allowed by the Super-Kamiokande event rates and the day- 
night spectra at la. It results in 0.97 < -ft'eorr ^ 1.03, which, as expected, does not differ 



much from the energy-constant value 1. Using Eq.(14) with K^^^ = 1.03 one obtains that 
the sterile oscillation is excluded with formal significance > 3.1 a. For the realistic evaluation 
significance is less (see Fig.l). 

If to add to the Super-Kamiokande data all other pre-SNO data (including the measured 
rates in the chlorine and gallium experiments), the status of sterile- neutrino oscillation be- 
comes further disfavoured. In this case the range of allowed values of Kcon is slightly shifted 
to lower values 0.93 < -ft'eorr ^ 1- With this range, the sterile-neutrino oscillation is excluded 
at significance better than 3.3 o", as it can be seen from Fig.l. 

For active neutrino-oscillation the analysis can be made in identical way. In this case the 
correction factors have to be calculated for all points in the parameter space (Arri^, tan^ 6*) 
allowed by all experimental data. In fact we used the range of calculated values of -f^eorr 



for the points of best fit solutions only (see Table 1) as 0.944 < -ft'eorr < 1.002. From 

Eq.(|l4|) we obtain formally rose = 0.244 ± 0.0695. From Fig.l we see that that observed 
fluxes for this range of correction factors are 3.3a away from rose = 0, which corresponds to 

hypothesis about absence of active-neutrino oscillation. 



^The analysis of sterile-neutrino oscillation, including the calculations of i^corr has been performed 
by M.C.Gonzalez-Garcia and C.Pena-Garay. They have also calculated the data of Table 1. 
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TABLE I. Correction factors Xcorr for the best fit solutions 



solution 




Am^ 


, eV^ 


K 

^ ^ corr 


LMA 


0.365 


3.65- 


10-5 


1.002 


SMA 


0.00061 


5.01- 


10-"^ 


0.944 


LOW 


0.708 


1.03- 


10-7 


0.993 




1.000 


5.46- 


10-11 


0.995 



In conclusion, we have derived the exact formula ([lID for calculation of ratio of the rate 
in Super-Kamiokande from i/^,- neutrinos to the total rate. This ratio, rose, characterises 
the oscillation of electron neutrino to another active neutrino. The formula is valid for arbi- 
trary thresholds in SNO and SK experiments, and arbitrary neutrino spectra. In particular, 
smallness of 5 Pee/ Pee is not required for validity of Eq.(|Tl]). We have shown that this ratio 
can be used to test the hypotheses of no-oscillations and oscillations to sterile and active 
neutrinos. 

For the no-oscillation case the ratio rose is given by Eq.(|13|) and this case is excluded at 
3.3(7 by comparison with theoretical no-oscillation value rose = 0. 

For a small energy- dependent component of survival probability, 5Pee/Peei the ratio is 
described with a good accuracy by the same Eq.(|TB|), which would exclude z/g Us oscillation 
at 3.3cr. This hypothesis is further tested evaluating numerically the correction factor -ft"corr 
for all points in the parameter space (Ar/i^, tan^ 9) in the pre-SNO allowed region for sterile- 
neutrino oscillations. It is found that sterile-neutrino oscillation is excluded at 3.3(T. 

For oscillation to active neutrinos our analysis is not as accurate as for sterile neutrinos. 
The correction factors are calculated only for the best-fit points in the parameter space. 
With these data the oscillations to active neutrino is confirmed at 3.3(7. 

With statistical error reduced by factor 2, the the significance of above results will reach 
3.5a. 
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